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Abstract. A new approach for obtaining the transformations of solutions 
of nonlinear ordinary differential equations representable as the compatibility 
condition of the overdetermined linear systems is proposed. The correspond- 
ing transformations of the solutions of the overdetermined linear systems are 
derived in the frameworks of the Darboux transformation technique. 



1. Introduction 

The nonlinear ordinary differential equations (ODE's) admitting a representation 
as the compatibility condition of the overdetermined system of linear equations 
(Lax pair) are intensively studied in modern mathematical physics. The most 
famous equations belonging to this class are the six Painleve equations Pi~Pvi 
|l| H, ^. It was found that these equations are closely connected with nonlinear 
partial differential equations (PDE's) integrable in the frameworks of the inverse 
scattering transformation (1ST) method [Q. Various approaches developed in the 
theory of nonlinear integrable PDE's were applied to study ODE's of such a class 
(e.g., see reviews ^ and references therein). 

It is well known that Pn-Pvi have the transformations that map the solutions of 
a given Painleve equation into solutions of the same equation but with different val- 
ues of the parameters of the equation. Such transformations for Pn were found in 
as a generalization of corresponding formulas of the rational solutions transforma- 
tions |8| . The transformations of solutions of Pn were also derived by means of the 
Backlund transformation technique . A representation of the Painleve equations 
as the systems of the first order equations was used to obtain the transformations of 
Pin in the cases of two different choices of values of the parameters jlo[ |ll| , Piv 
and of Py with special parameter values [pl]| . In similar manner the transforma- 
tions for Pn and Piv were obtained in Moreover, an equivalence of nonlinear 
ODE to a system of equations was used to construct the transformations of solu- 
tions of the higher order nonlinear ODE's appearing as the self-similar reduction in 
the KdV hierarchy ijlj, . An unified approach to investigate the transformation 
properties of Painleve equations, which utilizes the singularity structure of their 
solutions, were developed in |]l5| . The known transformations and new ones for 
Py, that differs from the equation considered in by a region of the parameters 
values, and for Pyi were derived. The truncation method was modified to find the 
transformations of nonlinear ODE's jl^. Its extension, which is based on consid- 
ering the transformations that preserve the locations of subsets of the singularities 
of solutions, was suggested The transformations of solutions of Pn and Pjy 
were, in particular, obtained using this approach. 
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The transformations of nonlinear ODE representable as the compatibihty con- 
dition of the Lax pair cause evidently the transformations of the Lax pair solu- 
tions. Studying the monodromy preserving deformations of matrix linear ODE's 
with regular singularities, Schlesinger [T^ constructed, in pure algebraic way, the 
transformations that keep the monodromy data of the solutions except for the ex- 
ponents of formal monodromy. These transformations, which are now referred to as 
the Schlesinger transformations, were generalized for linear equations with irregular 
singular points by Jimbo and Miwa The monodromy preserving condition is 
equivalent to having a set of additional linear ODE's on solution of the isomon- 
odromic equation which establish a hierarchy of the Lax pairs. However, 

the authors mentioned seem to apply no the Schlesinger transformations for the 
derivation of transformations for nonlinear ODE's (such as Painleve equations; see 
|3|, p. 437) that admit the compatibility condition representation. The transfor- 
mations of Pii-Pvi, including new ones, were obtained by means of the Schlesinger 
transformations of associated Lax pairs in |Q . 

In the present paper we show that the Darboux transformation (DT) technique 
|2l| can be used to construct the transformations of solutions of nonlinear ODE's 
and associated Lax pairs. This technique is suitable for generating the infinite hier- 
archies of solutions satisfying the reduction constraints imposed on the coefficients 
of the Lax pairs. To extend the DT technique for obtaining the transformations 
of ODE's we suggest that the points, in which the DT possesses a kernel, coincide 
with the singular points of the Lax pairs. The cases of Pn and certain nonlinear 
ODE of the second order are examined by way of illustration. 

The paper is organized as follows. The Lax pairs of the ODE's under consider- 
ation are written in Sec. 2. We present the explicit formulas of the DT's in Sec. 3 
and discuss the properties of them. The asymptotic expansions at the singular 
point neighborhood of the Lax pairs solutions are given in Sec. 4. These expansions 
are used there to build the transformations of the ODE's on the basis of the DT 
formulas. 

2. Lax Pairs and Nonlinear ODE's 
Let us consider direct Lax pair 

for matrix function = ^'(a;. A) and dual Lax pair 



H. . -SQ 

for matrix S = A). The compatibility condition of the direct pair ^'j;a — '^\x 
(or, for dual one, 'E.^x = Sax) leads to the condition on the matrix coefficients P 
and Q : 

Px-Qx + [P, Q] = 0. (3) 

The equation Pn 

Uxx = 2w'^ + xu — a (4) 
arises from the compatibility condition of Lax pairs (0j2|) with coefficients 

P = — A(T3 + U(7l, (5) 
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Q = 4A2(T3 - 4Aii(Ti + (73 {2u^(j-Y - 2v?E - xE) + X'^aai, (6) 

where Ui (i — 1, 2, 3) are the Pauh spin matrices, is 2 * 2 unit matrix. 

In the sequel the DT technique will be applied for obtaining the transformations 
of the solutions of Pn and following ODE 

x^uuxx — (xUx)'^ — xuux + cP'u^ — 2ju^ — 2(2/3 + l)xu ~ Ax^ (7) 

(a, (3 and 7 are constants), for which the transformations were not written explicitl 
in the literature. This equation is connected with special cases of Pm and Py [0, ^ 
Let us consider direct and dual Lax pairs (1,2) with coefficients 



X 



Q = --3 + A-f ^ \)+\-'{ V (9) 

\ V2 -p J \ W2 -q J 

v(cx — z) Q; — z (X 

where vi = uzw, V2 — , wi — zw, W2 — , <1 = z — —. The com- 



patibility condition (||) in this case yields 



2 



= 2 + ^z; - 2z^„2^ 

Z'x 2^ z{^z Q;), 
Wx _ uz + v{z - a) 
w ~ X 

The first three equations are reduced to Eq.(|^). 

3. Darboux Transformation Technique 

The first equations of Lax pairs (|l|,||) with coefficient P given by Eq.(|) or Eq.(|) 
are well known in the theory of 1ST as direct and dual Zakharov-Shabat (ZS) 
spectral problems 

= (AP(i) +P(°))*, (10) 

5, = -S(AF(i)+FW). (11) 
Potential P*-"-* of these problems has following general form: 

p(0) ^[0 ui 

{U2 

The DT technique allows one to produce new solutions of the ZS problems and 
corresponding matrix potential, starting from the initial ones. Let tp — ((^1,(^2)^ 
be a vector solution of direct ZS problem (|l^) with X = fi. The first elementary DT 
(EDT) of direct problem {*,S,p(")} ^ {^,E,P'-°^ is defined as given |2|, |2|: 

Ul(p2 

^=^1 (^^^^ I*. p, ^ I, (12) 
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p(0) 



2^ 



2 



2^1 




flUi 



(13) 



(it is assumed hereafter that cr is a scalar function of A). The direct and dual ZS 
problems are covariant with respect to this EDT: matrix functions ^ and S are the 
solutions of Eqs.( p^ , pl] ) with potential P^^K The second EDT of direct problem 
has a form: 



1 















2(^2 / 






/ 









ulipi 






2¥'2 



^ / A — /i — 



U2 

2 



-2^ 

^2 







(14) 



(15) 



In similar manner one can define the EDT's of dual problem, which depend 
explicitly on a solution of the dual ZS problem (pi]). Successive carrying out of an 
EDT of direct problem and proper one of dual problem leads to so-called binary 
DT (BDT) H ||: 



^crl E- ^ -R I S 



a \ fi — A 



R 



(16) 



p(o) ^ p(o) ^ _ ^) Tpd)^^- 



(17) 



where matrix R is the projector: Rij = fiXj / {XiV)i 'V — is^iiV^)^ is the vector 
solution of Eq.(|o|) with X = fi, x = (xii X2) is the vector solution of Eq.(pT[) with 

It is seen that the potential in the Lax pairs for Pn is the symmetric matrix. The 
DT technique is convenient for keeping the reduction constraints imposed on the 
coefficients of the spectral problems |^ . We can obtain the transformed potentials 
to satisfy this reduction, performing the iterations of BDT [|2^, [2^ . 

The DT's presented in this section allow us to build infinite hierarchy of solutions 
of nonlinear PDE's integrable in the frameworks of the 1ST method. It is supposed 
in so doing that vectors ip and x the solutions of Lax pairs associated with 
given PDE. Unfortunately, this approach is unfit for obtaining the transformations, 
which generate infinite hierarchy of solutions of nonlinear ODE's representable as 
the compatibility condition of Lax pairs. Vectors (p and x cannot be regarded as 
the solutions of systems (|[||) , since any transformation of the second equations of 
these Lax pairs ought to have no the kernel. The basic idea of the generalization 
of the DT technique for ODE's is to put the points, in which the DT of the first 
equations of corresponding Lax pairs has the kernel, into the singular points of the 
second equations of Lax pairs. For this aim we will assume in the next section that 
if and X s-re the solutions of systems (0,^ in the points X = fi and X = v. Then 
we will consider the limits in the formulas of DT's, tending the points /i and ly, in 
which the kernel exists, to singular point A = 0. 
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4. Transformations of ODE's 

To fulfill the procedure suggested in the previous section we need the explicit 
expressions of the asymptotic expansions of solutions of Lax pairs (|l|,^ with coef- 
ficients given by Eqs.(||,||) or Eqs.(||^) at singular point neighborhood. 

If a is unequal to half-integer, the asymptotic expansion of the Lax pairs solu- 
tions of Pn at neighborhood of point A = have the form 

^ ^>„{E + \p + o{X))A, (18) 
E ^ A-\E - Xp + o{\))%\ (19) 

Here we use notations 



-1 



The expansions in series of solutions of Lax pairs contain the logarithmic terms in 
the case of half-integer a. 

The asymptotic expansions of the solutions of Lax pairs of Eq. (|^) at neighbor- 
hood of point A = are given by Eqs.(^ 19), in which the coefhcients are defined 
as follows: 



\ WZ W ' 



(e( 

V eV2A a 



2\ a 1 



(21) 



where functions Ci and C2 solve equations: Ci^x = x ^ """^i' "^2,2: = ~^uc2. 

These asymptotic expansions allow us to apply formulas of the DT's for deriving 
the transformations of ODE's. Having the reduction constraint on potential 
we start from BDT (^6|jl^) in the case of Pn. To satisfy condition Sp Q = under 
performing the transformation we put a ~ 1. Substituting 

.=»u,(::), (22) 

X = (!>1, ti) H|a=„ (23) 

(oi, 02, 61, 62 are constants) and considering the limits p —t Q and — )■ in the 
BDT formulas, taking into account Eqs.( [T^JT9| , |20| ), we obtain two transformations 
of Lax pairs solutions and corresponding well-known transformations of Pn (e.g., 
see [H): 

* = ^' + ln:,(Al)(a3+^c^2)*/(2A) ( * = * + ln^(A2) ((73 - i(T2)^'/(2A) 
S = S-ln,(Ai)((73 + i(T2)S/(2A) I S = S - ln,(A2) (c73 - «fT2)S/(2A) 
u ^ u-\- h\x Ai ' I u = u — h\x A2 

a — \ — a y a = —1 — a 

The formulas of transformations coincide with presented above or follow from 
them as limits in the case of half-integer a. 

The transformations of Eq.(|^) are constructed similarly to that of Pn. Substi- 
tuting the asymptotic expansions (p^pl|) in formulas of EDT's (|l^,|l5) accordingly 
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to Eq.(J2^) and considering limit /i — > 0, we obtain four transformations of solutions 
and coefficients of Lax pairs: 






q- 


V1V2 
2x 


^2 ^2 ' 




At 




\ V2 - 

X 



P = P - 7^, vi = ~2xsf, V2 = Af--^, wi = vi + 2Pef - V2st^ 

W2 = 

1 + 1 ± a + P (7-1 + _)_ /J ^2 + 

where = ^ = , A^ = wi + -vi + —e^ , A^ = W2 + -^2 - tt^ 

ef 2wi X 2x X 2x 

(To keep condition Sp Q = we assume a = 

in the formulas of EDT's.) 
These transformations lead to the following ones of solutions and parameters of 

Eq.(0) {u, a, /3, 7} ^ {u, a, ~P, 7} {i = 1, 4): 

^ . 2x{xu^ + au^ - 2a; - 2{[3 + , „ 1 , a 

Ti- u = — ■ 2 ^"1 ^ — a = a, = 7 = 7- ; 

^ u(Q;(a;Ua; — au^ — 2a;) + 271*) ' ' 2' 2' 



^3 • 



2x(a;u2, + au^ + 2a; + 2j3u) _ n n ^ ~ ^ 

a = a, p = p — 7 



u{a{xux + au'^ + 2x) — 2ju)^ ' 2' 2' 

^ M(Q;(a;Ma; ~ au^ + 2x) + 27U) ' ' 2' 2 

The EDT's of dual pair give the same transformations. The transformations 
obtained commute and transformations T4 and T3 are inverse to Ti and T2 re- 
spectively. One can note that Eq.(^ is invariant under changing a — > —a. This 
symmetry connects transformations Ti and T2, T3 and T^. It should be stressed 
that a form of these transformations differs from one considered in |l5| . 

If 7 = ±a/3 or 7 = ±a{P + 1), particular solutions of Eq.(^ satisfy the Riccati 
equations. In these cases we obtain, performing the transformations of Eq.(^, 
which keep this condition, two transformations 

-ijj = y/x - ^ V\ { ^ = \fx + 



X 

p = p+l [ ~^ = p-\ 

of the linearization of mentioned Riccati equations: — ( K ^ — -] ^ = 0. 

\ X x^ ) 

So, the modification of the DT technique presented here allows us to obtain the 
transformations of solutions of the nonlinear ODE's admitting the compatibility 
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condition representation. At the same time we derive the Schlesinger transforma- 
tions of associated Lax pairs. This approach can be apphcd for finding the trans- 
formations of nonhnear systems of ODE's and their Lax pairs that have different 
types of the singularities and the reduction constraints of the coefficients. 
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